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Abstract 



A renormalization scheme is suggested where QCD input parameters - quark 

•/^ ' mass and couphng constant - are expressed in terms of gauge invariant and 

infrared stable quantities. For the renormalization of coupling constant the 

quark anomalous electromagnetic moment is used; the latter is calculated in a 

two loop approximation. Examination of the renormalized S matrix indicates 



P\J I confinement phenomenon already in the framework of perturbation theory. 
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Usually, in the charge and the mass renormalization procedure the so-called MOM scheme 
is used, where the input parameters of the theory, appearing in action, are expressed in terms 
of on mass-shell Green's functions (Itzykson and Zuber, 1980; Ramond, 1989). Another widely 
accepted scheme, MS, deals only with the divergent parts of the Green's functions, not appealing 
to any condition on momenta. Evidently, due to the renormalization invariance, any scheme is 
admissible. 

In gauge theories like quantum electrodynamics (QED) and quantum chromo dynamics 
(QCD) the input parameters are considered to be gauge invariant and infrared stable (i.e. 
not containing infrared divergences, generated by a massless gauge field). Let us for brevity 
call any quantity with these properties a physical quantity. Obviously, the observables, being 
measurable, should be the physical quantities but the converse may be not true (as an example, 
consider any function of field strength F^y in QED). Since Green's functions are not infrared 
stable and depend on a gauge (in general, their renormalization factors are nonlocal quantities 
(Basseto et al, 1987)), the interpretation of the results obtained from the schemes mentioned 
above may be obscured. Renormalized quantities may share undesirable properties, originated 
from the dynamics of the nonphysical degrees of freedom and the extraction of the physical 
information may be complicated. 

Therefore, from our point of view, in gauge theories the most illuminating would be a scheme 
operating only with the physical degrees of freedom, thus allowing to avoid complications 
mentioned above. The prescription is as follows: calculate in regularized theory as many 
physical quantities as the input parameters are and express the latter in terms of physical 
quantities. So, in QED as well as in QCD, containig only two input parameters - coupling 
constant and the fermion mass (for illustrative purposes we consider here only one flavor) - we 
need the expressions for two physical quantities: 

Sa = Sa{go, mo] n);IJ/3 = Ep{go, m^] n), (1) 

where Qq and rriQ are input parameters, the dependence of E on momenta, beinge irrelevant 
here, is omitted, n is a. space-time dimension and throughout this paper we use dimensional 
regularization. In QED and QCD Ea and Z"^ have no limit at n = 4 - regularization can 
not be removed in (1). Now resolve go and mo in terms of Ea, Ep and n and substitute 
into the expression for any other physical quantity E^. In renormalizable theories, after all 
the calculations are performed, E^ becomes finite in terms of E^, Ep - regularization can be 
removed: 

\mi^E^f{go,mo]n) = hm E^{go{Ea, Ep;n),mo{Eo„ Ep;n);n) = 

lim E^{Ea, Eb; n) = a^ < oo (2) 

This scheme, describing the renormalization procedure as the expression of physical quantities 
in terms of physical quantities,^ allows to avoid any significant difficulties and seems most 
transparent from the physical point of view. 



^A similar procedure (in a different context) was discussed already in the early days of QED (Dyson, 
1949). 



The aim of this paper is to develop such a scheme for QCD. 

It is clear that we need two physical quantities. In (Tarrach, 1981) it was demonstrated 
that in the covariant gauge the solution of the equation 

G~^{p, mo, go, ^; n)^i„(p) = (3) 

can be expressed in the framework of perurbation theory as 

m = mo + gl6mi{mo; n) + gQ6m2{mo; n) + ..., (4) 

where 5mi and 5m2 do not depend on a gauge parameter ^ and are infrared stable. In (3), G is 
a quark propagator and ^in{p) is the solution of the Dirac equation {■j^Pfj_ — m)'^in{p) = 0. The 
investigation in an axial gauge (Japaridze et al, 1991) confirms the gauge invariance and the 
infrared stability of m. In (Japaridze et al, 1991) it was shown that up to order gQ the quark 
propagator has a simple pole at p"^ = rrP'. Therefore, the quark pole mass can be considered 
as on of the Z"s in relations (1) and to complete the scheme we have to point out the another 
physical quantity in QCD. 

We propose the quark anomalous electromagnetic moment, defined, as usual, from the 
amplitude of the quark elastic scattering on an external electromagnetic field: 

(p|j^|p + A;)A^(A;) = 'l',„(p)('7M^i(A;2) + ^[7^,7,]rF2(A;2)U,„(p + A;)A^(A;) (5) 

The anomalous electromagnetic moment x is defined as -^2(0). 

We calculate x up to order g^ regularizing all the divergences (ultraviolet and infrared) in 
terms of space-time dimension n: 

X = 6-0 Xi("^o;n) +5-0X2 ("^o;n), (6) 

where Xi is the i-loop contribution. The gauge dependent terms cancel in Xi'-, Xi is infrared 
stable and the infrared divergence appears in X2- To obtain the expression of order g^, containig 
only 5o-associated divergences, we must use the relation (see (4)) 

mo = m - gl5mi{mQ] n) + 0(5-0) = ^ - gl5mi{m; n) + 0(5o) 

in (6), i.e. reexpand the loop expressions: 

X = glxiijn - gl5mi{m; n);n)+ 50X2 ("^; n) = 

gQXi[m;n)+ gQ\X2[m]n)-dmi — I (7) 

This reexpansion generates the infrared divergent tern dxi/dm which cancels the infrared 
divergent term in X2- 

So, X is gauge invariante and infrared stable. Omitting the intermediary calculations, we 
quote the final result: 



Cf5V. g^\llCA-2nf( 1 , ^' 1 ^ ^, 2 

+ ln—-\+^ + 0{gl-n-A) 



x = -^ 1 



Stt V 87r2 



^ 



n — 4 Attv"^ j 



(8) 



where Ca = N, Cp = {N"^ — 1)/2N are the invariants of SU{N) group, g"^ = fl'o'^""^ is the 
dimensionless couphng constant, u is the mass scale, appearing in the framework of dimensional 
regularization (Itzykson and Zuber, 1980; Ramond, 1989), n/ is the number of flavors, m is the 
pole mass of the scattered quark, the external momenta obey p^ = {p + k)"^ = m^ and for the 
finite part $ see the Appendix. 

Thus from the point of view of renormalization procedure QED and QCD do not differ - in 
both of theories the input parameters can be expressed in terms of physical quantities. Surely, 
any scheme can be used but our goal was to demonstrate that it is possible to renormalize 
QCD coupling constant in terms of physical degrees of freedom. The scheme is described by 
the relations: 

mo = mo(m, x; n), go = go{m, x; n) (9) 

To use (9) we need the numerical values of x and m. The gauge invariance and the infrared 
stability of these quantities does not mean necessarily that they are directly measurable, but 
these properties guarantee that the numerical values of x a-^^d m can be extracted from the 
experimental data. Of course, the numerical values of x a-^id m depend not only on mo and 
go but also on all other input parameters of, say, the Standard Model, but for the considered 
problem it is enough to analyze only QCD corrections. 

The scheme (9) may be not suitable from the point of view of numerical convergence - the 
rate of convergence depends on the numerical values of x and m. To improve the convergence, 
one has to introduce the so-called effective parameters (Itzykson and Zuber, 1980; Ramond, 
1989) gji and mn, i.e. to move to another scheme. It has to be pointed out that some statements 
formulated in terms of effective parameters (say, the increasing of the of the QCD effective 
coupling constant in the infrared region, sometimes interpretted as a physical effect of increasing 
the force between quarks at large distances) are scheme dependent and are not valid in another 
scheme. In other words, since the effective parameters are chosen arbitrarily, they can not affect 
the physical results. 

To see this, let us illustrate how the renormalization group equation and the renormalization 
scheme arise in a quantum field theory. It is transparent in the framework of dimensional 
regularization, where we have two parameters mo and g^ = (^qZ/""^. The mass scale z/ defines the 
dimension of go, providing the dimensionless action. Parameters g and u are not independent: 

g'o = g\i^iH"' = 9^^2)4-'', (10) 

i.e. 

^ + ^^/M = (11) 

The renormalization group equation (11) can be presented in a familiar form by introducing 
the effective parameter gniy) by means of the relation 

g^ = g\g\{v\v) (12) 

leading to 

.^=«..) (13) 



where 



(3 = lim 
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The renormahzation scheme is specified by the relation (12); then from (14) we obtain the 
appropriate /3-function. For example, the choice 
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(16) 



The schemes are labelled by particular choices of q (e.g MS is obtained if we choose q = 
aijin — 4)). The introduction of vriBiy) is based on the same argumentation. 

So, the behavior of the effective charge defined through any particular scheme (e.g. leading 
to asymptotic freedom), being scheme dependent, may not lead to any valuable results - the 
behavior and numerical values of effective parameters depend on our choice and are not defined 
from the theory alone. 

Thus, the advantage of scheme (9), besides gauge invariant and infrared stability is that 
it does not operates with the effective parameters, allowing us to avoid conclusions which are 
insignificant from the physical point of view. 

As it becomes evident, the behavior of g^, and vhb, does not lead to a scheme- independent 
statements, e.g. the absence of quarks and gluons in the asymptotic states. On the other hand, 
from the renormalizability ofQCD it follows that in the expansion of any physical quantity a^ 
(say, Wilson Loop) 



o^ 
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the coefficients a^^j are gauge invariant, infrared stable and contain no ultraviolet divergences 
at n = 4, i.e. a^^ j are finite. 

So, the following question arises: are the quark and gluon degrees of freedom observable, 
i.e. do quarks and gluons appear in asymptotic states? 

The answer may be obtained from the examination of the scattering matrix. We suggest 
the following criterion: if at least one S'-matrix element built up in terms of fields is finite, the 
appropriate quanta appear in asymptotic states as a particles. 

In QED it is well known that the electron elastic sacttering amplitude is infrared divergent 
and taking into account the emission of photons leads to the cancellation of these infrared 
singularities - only the inclusive cross sections are finite (Itzykson and Zuber, 1980; Ramond, 
1989). That is why we can say that from the Dirac-Maxwell equations it follows the existence 
of electrons and photons as an observable particles. Of course, we a priori know that they exist 
and the S'-matrix analysis is in accordance with the experimental data. 

We consider the S'-matrix element of scattering of quark on an external electromagnetic 
field. Note first that the existence of a physical quantity x — -^2(0) (see (5)) does not mean at 



all that the amplitude of the elastic sacttering is finite - in a full analogy with QED, the infrared 
divergences remain in the elastic amplitude after the charge and the mass renormalization is 
performed. Let us consider the gluon emission, assuming that the inclusive cross section might 
be finite. According to the LSZ reduction technique (Itzykson and Zuber, 1980; Ramond, 1989), 
for a gluon with the momentum q and the polarization e^(g) in an asymptotic state we have 

Z3 Z3 y 

1/2 
where V^p is the gluon propagator and Z^ is the gluon wave function renormalization factor. 

In order g^ the contribution of the gauge field in the residue Z3 is (Itzykson and Zuber, 1980; 

Ramond, 1989) 

^3 (O - ^C'^^^^TvV^ ^^37 T^;;^) (^) ' ^^^) 

where F is the Euler's function (Bateman and Erdelyi, 1973). The factor Z^ contributes to the 
QCD coupling constant renormalization. To obtain the expression for the amplitude we have to 
perform all the calculations before the regularization is removed: renormalize mo and go using 
(9), use (as we did for the amplitude (5)) the conditions p^ = (p + k)"^ = rn?, q^ = and then 
set n = 4. The condition g^ = means that Z^{0) can be equated to zero. This is analogous to 
the procedure of vanishing of integrals corresponding to a tadpole-type diagrams in a massless 
theory- because of the analyticity of theory in n we can always find the region where the result 
is well-defined and then continue analytically in the desired region of n (Itzykson and Zuber, 
1980; Ramond, 1989). So 

Z.,{q^) = Z^{q^) + Z/^'^™°"(g2) ^ z|™"(0) (20) 

at q^ — > 0. This means that only part of the go renormalization constant, namely Cyi/2 — 2?2//3 
is restored (compare to (IIC^^ — 2?2j)/3 in expression (8)). In other words, ^3(0) does not 
renormalize go. 

Therefore, if we consider the gluon emission in order to cancel the infrared divergences of 
the elastic scattering amplitude, the ultraviolet divergence appears and the cross section in 
order g^ contains an unavoidable infinity. The S'-matrix is the same in any scheme but the 
result is transparent while using the scheme (9), manipulating only with the physical degrees 
of freedom and not using the effective parameters. 

Though the discussuin above is not a proof, it can be considered as an indication on a 
confinement phenomenon already in the framework of perturbation theory. In other words, 
QCD might be an example of a field theory where fields do not necessarily refer to the physical 
particles 0. 



^ Long time ago Schwinger (Schwinger, 1962), starting from the idea that fields are more fundamental 
entities than particles, realized this conjencture in a two-dimensional electrodynamics, where the 
particles, corresponding to fermionic degrees of freedom, are absent in asymptotic states. 



The next step would be the consideration of the scattering of colorless bound states. At the 
present time we have no definite result for this problem. 

To conclude, in QCD it is possible to define renormalization procedure operating only in 
the space of physical degrees of freedom. Though the relations (9) guarantee that the physical 
quantities built up in terms of quark and gluon fields are finite, this is not enough for the 
existence of these field quanta as physical particles. The examination of the quark scattering 
amplitude indicates that quark and gluon field quanta do not appear in asymptotic states. 
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APPENDIX A: 



The finite part $ is 
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9 



- 8Ye + 27r^/n 2 - 3C(3) + 



12 



2^/7e ^ 49rv _ 22 ^ 



267r2 



+ 47I + CA 



'I3l7r2 I6757I 177_ 



959 



n 



1/2 



A 
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|/n A, -y)F.3 + \fI^ + (3 - 2/n A,)F,4 - 2FU 



(Al) 



where tti is the mass of the scattered quark, Aj = mf/m'^, rrii 7^ m, 
•je — 0.5771 is the Euler constant and 



00 1 

c(3) = E -^ 

j=lJ 



(A2) 



is the Riemann zeta function. 

The Fi/s are the generahzed hypergeometric functions (Bateman and Erdelyi, 1973) of Aj 
and the parameters of Fij are as follows: 



1 n+2 n— 3 n—2 



(A3) 
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and 



dF- 
jpi — "-^ »j 

'^~ dn 



(A7) 



The functions Fij and F^'- are considered at n = 4. 
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